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An intensive study for both the weak coupling and strong coupling limits of the ground state 
properties of this classic system is presented. Detailed results for specific values of finite N are given 
and from them results for general N are determined. We focus on the density matrix and concomi- 
tantly its Fourier transform, the occupation numbers, along with the pair correlation function and 
concomitantly its Fourier transform, the structure factor. These are the signature quantities of the 
Bose gas. One specific result is that for weak coupling a rational polynomial structure holds despite 
the transcendental nature of the Bethe equations. All these new results are predicated on the Bethe 
ansatz and are built upon the seminal works of the past. 

I. INTRODUCTION 

The one dimensional delta function Bose gas is a classic in the field of exactly solvable integrable systems. Following 
on the realization that in the infinite coupling limit, the impenetrable limit, the system had many of the properties of 
the free Fermi gas in their seminal work, Lieb and Liniger |2j solved the model exactly. They derived the Bethe 
ansatz and Bethe equations and went on to solve for the excitations in the thermodynamic limit 0, • 

In this paper, we present an intensive study of the finite system, in the weak coupling and the strong coupling limits. 
Extensive analytical solutions are given for finite values of N, the particle number, and from them, the analytical 
solutions for general N are determined. Except for a couple of papers that give the excitations for a system of three 
particles Q and larger number of particles [a@ this is the first intense study of this finite N body system for finite 
coupling. We concentrate on the principle quantities that are the signatures of the Bose gas. 

After a preliminary Section ^ and Appendix [S] introducing the Bethe ansatz and Bethe equations and their 
properties, we give extensive analytical solutions, in Sections IIIII and llVI with the aide of Appendices iBl and IO for 
the density matrix and concomitantly the occupation numbers. In Sections and IVII we do likewise for the pair 
correlation function and concomitantly the structure factor. As such we build upon the seminal work of Lenard 0, 
the very important work of Jimbo and Miwa 8], the Leningrad group @, and our recent works 
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conclude the paper in Section IV 111 with some further comments following upon all these results. 

The occupation numbers and structure factors for the modes are the experimentally realizable signatures of the 
Bose gas. The spur to the recent, over the past decade, revival of active interest in this system is due to Olshanii 
|l3j , who pointed out how this system could be realized in nature and to the current experimental activity seeking to 
realize it |l4l Il5l Il6l Il7l Il8l Il9| . We refer to the introduction in our earlier paper for a discussion of the relevant 
physical parameters involved. 

As all the new analytical solutions in this work are base d up on the Bethe ansatz, it is appropriate to observe that 
this year celebrates the 75 th anniversary of Bethe's paper 20] in which he introduced the now famous ansatz in his 
study of the Heisenberg spin-chain. Ever since, it has been a golden key in unlocking the solutions to many exactly 
solvable integrable systems. Baxter [2l| gave a concise review of these in his tribute to Yang, who has recently written 
[22| on the Bethe ansatz to honor Bethe. 

II. CONSTRUCTING THE WAVEFUNCTION 

The purpose of this section is to present the Bethe ansatz wavefunction and concomitantly the Bethe equations for 
the one-dimensional delta function Bose gas in periodic boundary conditions, as derived in the seminal paper by Lieb 
and Liniger [2|. The Schrddinger equation for this system with h = 1, 2m = 1 is 



N ri 2 

^ dx'j 



i=i 



2c S(xi 

l<i<j<N 



Xj) ip N (x 1 ,x 2l ■ ■ - ,x N ) = Eip N (xi,x 2 , ■ ■ ■ ,xn) 



(2.1) 



where c controls the strength of the ^-function. Throughout this paper we consider only c > 0, the repulsive case. 
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Here we display the solution for periodic boundary conditions 



P&Sn 



(2.2) 



where Sn is the symmetric group on N symbols, and hence the wavefunction is a sum over N\ quantities. The 
function e(p) is the signature of the permutation p, and the function a{p) is given by 



N 



»(p)=n 



1 + -( k PW -K(i)) 



(2.3) 



Note that a(p) may be given in a number of forms (any change absorbed into the normalisation), we choose <|2.3[) 
from Jimbo and Miwa [8j . Throughout this paper we use the normalisation given by 



N 2 



dxi... dx N -i \ipff(0,xi, x N -i)\' 
Rn—i 



(2.4) 



where Rn-i is the domain of integration specified by 



R 



N-l 



< x\ < . . . < xjv_i < L. 



(2.5) 



In this paper we are concerned only with the ground state of the system, where J^Li ki = 0, which leads to 



+i-i 



Vt = 1, 



, N. 



(2.6) 



The ki are ordered such that fc^ > fc^r-i > ... > ki > k\. The N (real) numbers ki are determined as the solution 
of the Bethe equations. These can be given in many forms 0,11,123]' we display that of 



A' 



(k j+ i - kj)L = 2tt + Y^ 



i.=i 



2 arctan 



kj k 



2 arctan ' 



j = l,...,N-l. 



(2.7) 



While these equations cannot be solved explicitly for kj as a function of c, they can be solved for both small and 
large cL expansions, using the method of quadrature. We list explicit small cL expansions in Appendix ^ and at 
this point we highlight that 



AN) 



24 



(cL) + 0{cLf 



(2. 



where hj* is the jth zero of the TVth Hermite polynomial. The leading term in (|2.8|l is found in reference to this 
problem by Gaudin |2^, for more detail refer to Szego |24|. The coefficient —1/24 in the next term of the expansion 
is new, and appears to be universal (see Appendix lAl . 

Here we have chosen to display the large cL expansion for kj in its general N form as given in 8] , for its particular 
utility in the following sections. 



(2J-JV-1)- 



1 - 2N 



cL 



AN 2 



cL 



+ \ -8N 3 + —N [2 j 2 + (N + l)(N — 2j)] vr ; 



(2.9) 



We also list some exact solutions to the Bethe equations in Table [I] 
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N 


cL 




2 




?r/2 


3 


tt/2 


7r(vTf-3)/4 


3 


TV 


V2tt 


3 


3tt/2 


3tt(v / 17 + 3)/4 



TABLE I: A sample of exact solutions to l|2.70 



It is now possible to explicitly construct the wavefunction iPm{x\ : . . . ,xn) using Ij2.2|l and We close here by 

exhibiting the unnormalised wavefunction ip2(xi, x 2 ) by way of example 



M*u*2) = ( 1 - — 1 e ik *^-^ - ( 1 + ^ ) e -^(— 0. ( 2.io) 



c / \ c 



which has normalisation 



sin2fc ? i 8sin 2 fc 9 i 4k 2 (2k 2 L + sm2k 2 L) 
M 2 =2L — + v 5 — . (2.11) 

fc 2 c c z 

In all that follows the weak coupling expansion corresponds to the dimensionless parameter cL < 1 and the strong 
coupling expansion corresponds to cL 3> 1. 



III. DENSITY MATRIX AND OCCUPATION NUMBERS 



The normalised density matrix pn(x, 0) is defined as 



N 1 



/ dxi . . . dxjv-iV'iv(0,xi,X2, 

■_ JR N - 1 , j (x) 



Pn(x,0) = — j-p 2^ I dxi . . .dx N -itp N (0,xi,x 2 , ■ ■ ■ ,x N -i)tp N (xi, . . . ,Xj,x,x j+ i, . . . ,xn-i) (3.1) 
3- 



where the overbar implies complex conjugation, the normalisation Af 2 is specihed by (|2.4(l . and the domain of inte- 
gration is specified by 



Rnj(x) : < x\ < . . . < Xj < x < < . . . < xn < L. (3-2) 

The density matrix is normalised such that Pat(0, 0) — po — N/L. Hence to compute the density matrix for A" 
particles, one must perform N lots of N — 1 dimensional integrals over (AM) 2 terms (e.g. 4 triple integrals over 576 
terms for N = 4). This is a computationally expensive task, and hence we were able to determine solutions only up 
to A^ = 4 using this method. 

The occupation numbers c n (N) are determined as a Fourier transform of the density matrix 



c n (N)= p N (x,0)e 2i7Tnx/L dx = / p N (x,0)cos(2irnx/L)dx. (3.3) 
Jo Jo 

They have the physical interpretation of being the expectation value of the number of particles in mode n. 

We note the normalisation property X)n=-oo c «(-^0 = we have confirmed this result for all occupation number 
formulas that follow. It is sometimes useful to discuss occupation number per particle; we introduce the notation 
c* (N) = Cn(N)/N. We shall now display some explicit solutions for N = 2, 3,4 and general N. 
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A. N = 2 



Within Section ITTll this subsection heralds the most complete set of results, with results for N = 3 and 4 becoming 
increasingly exiguous as the intricacies of the equations develop. For example, it is only possible to display the 
complete density matrix for N = 2, as for N = 3 already the equation would take many pages to display. We have 
utilised for N = 2 (|2~T|) to obtain 



k 2 L 

• = 2k 2 tan [ — — 



hence producing a concise form of the density matrix using l|2.2|l . (|2.3|l and (|3.1|l 

2 k 2 x cos(k 2 (L — x)) + k 2 (L — x) cos(k 2 x) + sin(fc2 (L — x)) + sin(fc22:) 
P2(X ' } = I fc 2 i + sm(k 2 L) 

with corresponding occupation numbers from (|3.3I) 

4(fc 2 L) 3 (1 - cos(k 2 L)) 



c„(2) 



' (4n 2 7r 2 - k 2 L 2 ) 2 {k 2 L + sin(fc 2 i)) 
The expansion of the density matrix for small cL is given using IIA1I) 



< k 2 < TT. 



< fc 2 < 7T 



(3.4) 



(3.5) 



(3.6) 



p 2 (*,0) = - 



i- f2( :,- f)2 (cL)' 



t 2 {ll-t) 2 {t 2 -Tit + 27T 2 ) 



24tt 4 v ' 360tt 6 

t 2 (7T - i) 2 (167T 4 - 247T 3 i + 277r 2 i 2 - 67rf 3 + 3i 4 ) 

40320tt 8 



(cLf 



(cL) 4 + 0{cLf 



(3.7) 



Note that we have introduced here t = ttx/L, henceforth we switch between t and x as appropriate. The corre- 
sponding occupation numbers for l|3.7l) are given by 



C„(2) = 2 ' 1 _ 720 + 6Ms( cI/ ) 3 ~ 1209600 ( cL ) 4 + °( cL )" 



when n = 



16 



4Tt 4 7T 4 -30n 2 7T 2 +45 



:!sin 



+45 (cL) 4 + 0(cL) 5 when n ^ 0. 



(3.8) 



Utilising (|2~3|) . we obtain a large cL expansion for the density matrix 



P2M) 



2 r (it -2t)cost + 2sint 8(ir - t)tsint f 1 



s |(,t 2 - :-!,t/ + 2/-)/ cos / ~ (V- -! (S;t/ - - sint] ^ 1 \ 2 + Q / 1 \ 3 j 



(3.9) 



Note that in the limit cL — > oo, we recover (20) from [lOj . The corresponding occupation numbers for (|3.9(l are 
given by 



c„(2) = 2 



32(12ti 2 + 1) / 1 



(4n 2 -l)V (4n 2 -l) 3 7r 2 \cL 
{ 3072n 2 (4n 2 + l) 64 (6n 2 



1 



(4n 2 - 1) 4 tt 2 



(4n 2 - l) z 



1 



O 



cL 



(3.10) 



which in the limit cL — > oo recovers (42) from [T^. We also display here the density matrix for the exact solution to 
(|2~7|) from Table|U when cL = tt, fc 2 L = tt/2 



P 2 (i,0) 



2 (-t + tt + 2) cos | + (t + 2) sin : 

L 7T + 2 



(3.11) 



5 



and the corresponding occupation numbers are given by 



c n (2) = 2 



16 



(16n 2 - If 7r(7r + 2) 



(3.12) 



B. N = 3 

The small cL expansion of the density matrix is given here, using ||OJ and (H}. 



1 12tt 4 {CL) 180^ (CL) 

t 2 (Tr - tf (— 101i 4 + 2027rt 3 - 1257r 2 t 2 + 247r 3 i + 54tt 4 ) 



■{cLf + 0(cLf 



(3.13) 



20160tt 8 

Note that the coefficient of the (cL) p term is a polynomial in t of order 2p, for p > 2. This structure is also repeated 
in lETTl . 

The corresponding occupation numbers for (|3.13|l are given by 



c n (3) = 3 



Ki.'i 



(cL) 4 + 0(cL) 5 



when n = 



^(cL) 2 - ^4#(cL) 3 + foT a, T V 1515 W 4 + °( cL ) 5 when » ^ °- 



(3.14) 



We give here the large cL expansion of the density matrix, using 13.1fl and 12.7fl 

, , 3 r2(7r-2i) 2 + 12(7r-2<)sin2£ + 4(2<-7r + 2)(2f-7r-2)cos2< + cos4i+15 

P ^ 0) = IX 6^ 

4sint [(7r- 2t) (l + 8?rf - 8t 2 ) cost - (tt - 2t) cos 3t + 8 (tt - £)tsini] / 1 



+ 



which in the limit cL — > oo recovers (21) of [13. 

The occupation numbers corresponding to (|3.15() are given by 



cL 



O 



Cn(3) =3< 



f (* + #) + (§ + £) (£) + o(£) s 
|-(! + ^)(^) + o(^) 2 



35 385 / 1 



1087r a 36 



when n = 
when |n| = 1 
when \n\ = 2 



2(3n 2 + l) _ 4(9« 6 -28n 4 -61n 2 +8) 



* ,,+ „- i.,-'.,,,-- i ! --' V ch) + ^ ( eh) when|n|>3. 



(3.16) 



which in the limit cL — > oo recovers (45) of jlO(. We also display here the density matrix for the exact solution to 
(j2~7jl from Table ID when cL = V2ir, k 3 L = tt 



rSt 2 - 8irt + 2(f>{ir - 2t) cos t - 4 cos 2t + 4 [tt (2y/2t + 5) - 2^2 (t 2 - 5)] sin * + <pn + 52 
A 48 + 3</>tt J (3 ' 17) 



LI 48 + 307 

where = 10\/2 + 3n. The corresponding occupation numbers for (|3.17|1 are given by 



_5_ , +\/2 _ 4(299+71 V2tt) 
27 ' it 27(16+0tt) 



, i 64x/2+54 7T 
Cn[o) — 6 < 8i7r(i6+<£7r) 



when n = 
when |n| = 1 
when \n\ > 2 



(3.18) 

4(-7T+32n 4 (2\/2+7r)-4ri 2 (l2\/2+7r)) 
3n 2 7r(4n 2 -l) a (16+0ir) 

The density matrix and occupation numbers for the other exact values listed in Table [I] for N = 3, although 
calculated, are too lengthy to display here. 
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C. TV = 4 

The intricacy of the density matrix at this value of N is already so great that we go on to calculate the occupation 
numbers without explicitly exhibiting it, and as such we give only Cq(4). We utilise (|3.3|) and (|A3|) to produce the 
result 

r 1 13 383 i 

D. General N 

It is interesting to observe that with the presence of the irrational numbers in the Bethe equations for small cL 
(Appendix^), that when the final results for the occupation numbers appear they contain purely rational numbers. 
Encouraged by this remarkable observation and other indications in the preceding subsections, we looked for a pattern 
in N for the coefficients in the small cL expansion for a general cq(N). Upon close examination of 13.6(1 . 13.14JI and 
(|3.19|) . and with some good fortune, we found the following polynomial structure for the small cL expansion of the 
n = occupation number for general N. 

01 ' 720 V ; 720 V 42 y 



N-ll /45(/V - l) 2 - 5(N - 1 



720 42 V 120 



{cLf + 0(cLf. (3.20) 



We have written (|3.20|l as given to emphasise the detailed structure of the coefficients in this expansion. Note that 
the coefficient of (N — l)(cL) p is a polynomial of order p — 2 in N — 1 for p > 2. Therefore to obtain say the coefficient 
of the (cL) 5 term we would need cq(N) for four specific values of N in their rational form. Unfortunately, it was too 
computationally expensive for us to obtain any N — 1 polynomial for any higher order than (cL) 4 . We conjecture 
that this is the pattern for all N > 2. 

The fact that the coefficient of (N — l)(cL) p appears to be a polynomial of degree p — 2 in N — 1 suggests that for 
any finite N there is always an interval cL E [0, Dn) such that the series is convergent, but with Dm — > as N — ► oo. 
To quantify this last point, note that in the thermodynamic limit the dimensionless parameter is c/po = r and the 
coefficient is proportional to N p , which suggests that the corresponding radius of convergence is proportional to 
1/N 2 . In particular this means that no information can be gleaned as to the functional form of Cq(N) as a function 
of r about r=o except that it is not analytic. 

IV. DENSITY MATRICES AND OCCUPATION NUMBERS FOR LARGE cL 

In the previous section, we gave large cL expansions for the density matrices and occupation numbers for N — 2 and 
N = 3. Again, it was numerically prohibitive to go beyond N = 4, furthermore as can be seen from the coefficients 
in these expansions (and as can be witnessed in the Tables ITlIIVIIII which we will refer to in what follows), that the 
numbers are highly irrational with no hope of finding an analogous pattern as we were fortunate enough to do in 

(£29. 

We therefore turn to a totally different mathematical strategy to obtain a large cL expansion for these quantities. 
In the impenetrable limit (cL — oo, arbitrary L), Lenard developed the theory for the density matrix for arbitrary 
N, and went on to show that for asymptotically large N that cq(N) ~ N 1 / 2 . 

In our recent work |l0j | we employed Lenard's theory to obtain the results for the occupation numbers for a range 
of finite iV, and from them determine the results for general N, which continue on to the asymptotically large N limit 

m 

To go beyond the impenetrable limit for these quantities, we turn to the very valuable work of Jimbo and Miwa 
Q . Building upon the work of Lenard , they developed an expansion for the density matrix in the large cL limit 
for general N in principle. We say in principle because while their expansion is superb in the form given, it is as 
numerically prohibitive to use as was the method we employed in the previous section. 

To resolve this difficulty we have recast their theory using mathematical techniques evidenced in our recent work 
[l2| and presented in great detail by Forrester [25| into a new form that is readily amenable to numerical calculation. 
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We will find in what follows that the specific results for N — 2 and N — 3 as in Section ITJIl are useful specific checks 
to the theory. 

In Subsection IIV Al along with Appendix FBI we present the full details of our derivation for the the density matrix. 
Following upon that, in Subsection IIV Bl we are now able to calculate the occupation numbers for a finite range of N 
values and from these results are able to determine the results for general N which again continue to asymptotically 
large N. 



A. Toeplitz Determinants 

The fact that \4>n\ 2 consists of (AH) 2 terms means any method based on term-by-term integration must necessarily 
be restricted to small N. To overcome this one must seek out structure in the form of ippf, and this structure must 
be used to reduce the computational expense required to compute pn{x,0). Certainly in the limit cL — > oo there is 
structure in l|2.2|l for then l|2.9fl gives kj — (2j — N — 1)tv/L while (|2.3|1 gives a(p) = 1, and so 



N 



,x N ) 



cL — >oc 



det[e?< 2k - N - 1 >j/ L ] j>k = lt .., tN . 



(4.1) 



(4.2) 



While the sum consists of N\ terms, the determinant can be computed in 0(N 3 ) arithmetic operations. Moreover 
the corresponding density matrix can also be expressed as a determinant (Lenard 7]) 



(4.3) 



where 




Xe -i*{N-l)x/L det 



3,fc=l,...,AT-l 



(4.4) 



Mj-k) = - 




2n\u/L g27rix/-£/Wg— 2-iriu/L ^^27riu(j — k) / L 



when j — k = — 1 
when j — k = 
when j — k = 1 



^ {j -k)({]-k) l >-i)-n l(j ~ k ) cos (0' ~ sin< ~ costsin((j - k)t)] when \j - k\ > 2 



(4.5) 



(4.6) 



and thus also can be computed in 0(A^ 3 ) operations. In [lOj this determinant formula was used to compute the 
density matrix and the corresponding ground state occupation numbers up to N = 7. It has been shown by Jimbo 
and Miwa || (see also Section IVTl below) that determinant structures persist if ipN is expanded in large cL, and that 
this implies special structures for the expansion of the density matrix. In particular, writing 



PN ( Xl 0) = p%\x,0) + (J^j p$(x,0) + (±J 



(4.7) 



it was shown in |8| that 



p { $ {x, 0) = -2 Po x—p { ° ] (x, 0) + F N (x) 



(4.8) 
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where (x, 0) is specified by (|4.3|) and 



F N (x) 



1 



A(x;-2) 
Introducing the kernel function 



IaP; A |a ( ,; A) + §-J J ; A ^(x; A) - A, * ; A ^ A(,; A) 



A=-2 



(4.9) 



1 W 



(2j-Af-l)(x-y)/L 

~ 3=1 

sin[A^7r(a; — y)/i] 



Lsin[7r(a; — y)/L] 



(4.10) 
(4.11) 



the quantity A]l ; A is defined in [8| as the Fredholm minor 



Ai *;A 



~ a ;+1 /■* 



1=0 

x det 



Knl{x,v) [ K NL(x,Uk)] k = Xt _ tl 



while A(x; A) is specified as the Fredholm determinant 



°° \l rx rx 

A(a;;A)=Y]— / dui . . . / dm det [K N , u k )}, k=1 
i=o 11 Ja Jo 



(4.12) 



(4.13) 



Neither l|4.12[l or (|4.13|) are suitable for computation. However for A x we have the determinant formula 14.4fl . and 
A too can be expressed as a determinant. This can be seen by developing Lenard's Q derivation of (|4.3|> . which 
proceeds by writing the summation of determinants of l|4.12[l and (|4.13(l as multidimensional integrals of a type which 
can be recognised as Toeplitz determinants (see Appendix IBl . In addition to (|4.4(1 - (|4.6(l . one deduces that 



A(x;A) =det [A (j - k)h fe= , ... 



N 



(4.14) 



where 



A (j - k) 



J^i^J + A J J duexp[2-!tiu(j - k)/L] 

{ l + ^f when j - k = 

I 2bF0=fcJ {e^^ k)xjL - 1) when j-k^0. 



(4.15) 



(4.16) 



With these determinant formulas Fn(x) is expressed in a computable form. This concludes the method necessary 
to construct the density matrix expanded in large cL. We now examine the occupation numbers. 

B. Occupation Numbers for Large cL 

The following notation for the occupation number c„ (N) expanded in large cL is given as 



c n (N) = c^(N) 



(4.17) 
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We list some exact values of c^ 0) (JV) and c4 1} (JV) for n = Q (JV = 2, . . . , 7) and n = 1, 2 (JV = 2, . . . , 6) in Tables ITTT1 
HV| and E| The c^JV) term can be expanded further as 

= I" (-2N X ^^ + F N (x)) cos ( ^) dx (4.18) 



^- + ^ W Jcos^— J 
= c^(N)+c^(N) (4.19) 

where Fn(x) may be computed using the expressions in Subsection IIV Al above. The Cn'^^N) term can be simplified 
using integration by parts to 

f L 

c£ A) (N) = -2JV 2 + 2Nc^(N) - AirnN / - sm(2irnx/L)p^(x)dx. (4.20) 

Jo L 

We list some exact values of ci M) (JV) and cl 1,2) (TV) for n = 0, 1, 2 with JV = 2, . . . , 6 in Tables EH E3J and IVlTTl 
We present all the numerical values that we computed for Cn°'(JV), c^ 1 ' i \N), Cn 1,2 ^(N), and (JV) obtained for 
n = 0, 1,2 in Appendix|0 In this table, the data is given to 6 significant figures for economy of presentation, while in 
point of fact, accuracy to 10 significant figures was needed for the analysis that now follows. We were able to achieve 
this numerical accuracy up to JV = 36 for the n = mode, and up to JV = 26 for the modes n = 1 and n = 2. Note 
that dn°'(N) for n = 0, 1,2 reproduces the results of [Icj . 

We now introduce a generalised version of the ansatz first introduced in [Tc| . expanded to 0(1/ cL) 



a n N \ ^i+iiiN ^ ( ^ 7„JV 



c„(JV) = A*,,,, ( 1 + -^-\ m + — + C oc ,„ I 1 + ^ ) (4.21) 



n „ / 1 ^ 2 



= (A»,nVJV + Coo, n ) + NiAeo^atnVN + A 00in /3 n VN]nN + C , 00in7n ) l-j + O l-j . (4.22) 

It is now pertinent to determine bounds on a n and /?„, and a value for j n . Consider the set of three linear equations 
describing the 1/cL term of l|4.22l) for three consecutive values of JV 

VN~T VJV - Tln(JV-l) l\ / ^oo,„a„\ / c£ (1) (JV-l) 

/JV VJVln(JV) 1 I I Aoo^Pn ] = I C ; (1) (JV) j (4.23) 

c ; (1) (jv+i) 

and the set of two linear equations describing the 1/cL term of l|4.22|l with = for two consecutive values of JV 



(4.24) 













H 




\ Coo , n J 





VJV VJV ln(JV) \ ( ^oo,„a„ \ = / C ; (1) (/V) 

VjvTi VjvTTMjv + i) ) { Aoo.nA, J I c ; (1) (jv + i) 



Values of Aoo^ and C^,™ from 01 are displayed in Table ITU The solution for a n , (3 n , 7„, of I4.23|) and the solution 
for a n , (3 n for (|4.24() for various values of JV establish bounds on a n and /3„, and a value for 7„. We establish numerical 
stability for these bounds by calculating them for N = 2 up to JV = 35. A fractional accuracy of > 10~ 10 is necessary, 
and hence the parameters are calculated at N — 35 for n = 0, and JV = 25 for n = 1 and n = 2. We present them in 
Table CU 

Note that f3 n is very close to 2 for n = 0, and suggestive of the value 2 for n = 1 and 2. We postulate that (3 n = 2 
for all n. We say more about this /3„ in Section IVTT1 
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A 

j^co ,n 


n 

woo ,n 


QL n 


a 
P 


1 


7« 





1.54273 


-0.5725 


0.1561 < qo < .1838 


1.998 < /So 


< 2.003 


0.1599 


1 


0.5143 


-0.5739 


-5.709 < ai < -5.067 


1.972 < fa 


< 2.094 


-1.109 


2 


0.3676 


-0.5775 


-8.350 <a 2 < -6.121 


1.887 < /3 2 


< 2.313 


-2.736 



TABLE II: Parameters for the ansatz I4.22H 



V. CORRELATION FUNCTIONS AND STRUCTURE FACTORS 



The definition of the two-point correlation function gN(x,0) is 



N-2 



9N 



( a; ' ) = X / dx 1 ...dx N - 2 \i[} N (Q,Xi,...,Xj,x,Xj +1 ,...,x N -2)\' 1 
■A/ Jr n ^{x) 



(5.1) 



where the normalisation A/" 2 is specified by (j2.4|l , and the domain of integration by (|3.2() . Then 1)5.1(1 has the property 



<7jv(x, 0)d.T = N — 1, 



(5.2) 



in keeping with the interpretation of gN{x, 0) as the density of particles at position x, given there is a particle at the 
origin. The structure factor is defined by 



S n (N) = 



1 1 



dxi . . . cLxn-i 



N-l 

E< 

3=0 



,2'iirxjn/ L 



\ip N (0,x u . . .jjcjv-i)!' 



(5.3) 



where xq = 0, which can be written in terms of gjv(a;, 0) to read 



S n (N) = 1 + / 3at(x, 0) cos(2irnx/ L)dx. 



(5.4) 



In view of l|5.2|l this gives S (N) = N. 



A. N = 2 



This subsection yields the most complete set of results within Section with the set of results for N — 3 and 
N = 4 becoming increasingly limited. For N — 4, we display only the small cL expansion of the correlation function 
and structure factor. 

It is possible to produce a simple, exact form of the correlation function for N — 2 by utilising (|3.4|) to obtain 



2fc 2 cos 2 [lfc 2 (£-2q : )] 
fc 2 L + sm(fc 2 L) 



which has corresponding structure factor 



2 when n = 

<Sn(2) = { 1 . k 2 2 L 2 sin(k 2 L) , , n (5-6) 
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The correlation function, expanded about small cL using IjAlfl . is given by 

, n . 1 r / x 2 x 1\ . _. / x* 2x 3 x 2 x 1 i r > 



L 



2x 6 2x 5 7x 4 x 3 7x 2 x 1 



45L 6 15L 5 36L 4 6i 3 90L 2 60L 945 
and hence the structure factor is given by 



{cLf + O(cL) 4 (5.7) 



. ( 2 when n = 

= I 1 - s£pr(<*) + T^(^) 2 ~ m4 'r^ a+60 (^) 3 + O(ci) 4 when n + 0. 
The large cL expansion of the correlation function is given using 



52(i,0) = -jr j 2sin 2 i + 8sin t[(7r-2i) cost- sin i] +0 | (5.9) 

taking the Fourier transform then yields the structure factor 

when n = 

5 n (2)= <( § + 3(i)+0(-U 2 whenH = l (5.10) 




4 (tt) +0(Jr) 2 when |n| > 2. 



2 1 \cL) ' V ci 

The correlation function for the special values cL = ir, k^L = ir/2, (see Table[Q| is 



and hence the structure factor is 



7r (sin t + 1) 

^o) = 4^W (5 - n) 



{2 when n = 

1 2 , n ( 5 - 12 ) 

1 - (4^-1)^+2) whenn^O. 



B. TV = 3 



Given the length of the full correlation function for N = 3, it is not possible to display the full result here, however 
we shall display the small and large cL expansions for both the correlation function and the structure factor. 
The small cL expansion of the correlation function is given, using l|A2|) 



93(x, 0) = - 



2 


H 


2 

X 


x 1 


L 


~1 2 + 


L ~ 6 




f X 6 


3x 5 


5x 4 


+ 


Ui 6 


5L 5 


' 8L 4 



X 

4L 3 



{ X X X X li-i 

\VU? ~ 6L 3 + 4L 2 ~ 6L + 40/ ( ° ' 



x 
40L 



1 

280 



(cL) 3 + O(cL) 4 



(5.13) 



and hence structure factor by 



5 B (3) 



r (cL) 



2n^7r -3 



(ci) 5 



nWlSn^Tr -180 f„r\3 
40n^<i 



when n = 
O(cL) 4 when n ± 0. 



(5.14) 
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The large cL expansion for the correlation function is given using l|2.9[l 



93 (i, 0) 



2 j 4(2 + cos2i)sin 2 i 8sint [-2 (tt - 2t) (2cosi + cos3i) + 3sini + 4sin3i] ( 1 



which yields the structure factor 



5„(3) = { 



1 , 32 
3 ' 9 \cL 
2 



_ + 38 /J_ 

3 ' 9 \cL 



) + o(^) 2 
) + o(^) 2 



when n = 
when |n| = 1 
when Inl = 2 



16(n -2) 



(n'-4)(n"-i) V^r) +°(^) when|n|>3. 



O 



cL, 



We now display the correlation function for the special value k^L = n, cL — \/2tt : (from Table HJ 

2?r [2 {Ay/2 - tt) cosi - 7rcos2i - 4 (2 + \72vr) sini + 2 (6 + \72tt) sin2i - 9n - 8y/2~\ 



<fe(*,0) = — 



and the corresponding structure factor by 



L (16 + 16\/27r + 97r 2 ) 



(5.15) 



(5.16) 



(5.17) 



3 

. i _ 32+16V27r-37r 2 
5 n (3) = ^ 48+48\/27r+277r 2 

16(2+V2tt) 
(4n 2 - 1) ( 16+16 </2tt+9it 2 ) 



when n = 
when |n| = 1 

when Inl > 2. 



(5.18) 



C. TV = 4 



Here we display the correlation function for the small cL expansion of N = 4 



1 + 



( 



x 2 x 1 
1? + L ~ 6 
7x 6 7x 5 47x 4 



(cL) 



4 3 1 

X X X 1 

W + 3L 3 " ~ 6L + 30 



(ci) 2 



2a; 3 19x 2 



1 



V18L 6 6L 5 36L 4 3L 3 
and its corresponding structure factor 



180L 2 30L 135 



(cLf + 0(cL) 4 



(5.19) 



5„(4) 



when n = 



^(ci) + ^_M( ci ) 



2 " V4 + 6 °gy" 525 (^) 3 + Q(^) 4 whenn^O. 



(5.20) 



D. General N 



Through close examination of the small cL expansions of the correlation function for N = 2 1)5.7)) . = 3 (|5.13|) 
and TV = 4 l|5.19[l . we find the following polynomial structure 



g N (x,0) = ^i{l-/(z)( c i) + 



-4+6 + (l2-9 / ^ + (720 - 216 



(cX) 2 



.^ + 23iv_ n /(T)3 

36 60 10 J J y ' 
144 720 120 J y ' 



36 40 + 5 I J[ ' 



/v- 



797V 



1 



6804 90720 1080 



(cL) 3 + 0(cL) 4 } 



(5.21) 
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where f(x) — x 2 / L 2 — x/L + 1/6. We conjecture that this structure continues for all JV. This has corresponding 
structure factor 



S n (N) = 



N 

(A'-l)(jVn 4 ir 4 +75A r n 2 7r 2 -180n 2 7r 2 +75A r2 - 1035^+1890) 



when 



(ci) 3 + O(cL) 4 when n 7^ 



(5.22) 



We have not proceeded further than JV = 3 for the large cL expansion for the same lack of utility as we encountered 
in the case of the density matrix in Section ITTT1 There is now, once again, a powerful way to proceed to develop such 
an expansion for general JV, and we present the derivation in the following section. The special case of JV = 2 Ij5.9|l 
and JV — 3 l|5.15|l provide useful checks. 

In closing this section, we note that in i|5.21|) . the leading order term for g^(0,0) is (JV — l)/£, which is the free 
Bose result. Deviations from this begin at order cL. 

VI. CORRELATION FUNCTIONS AND STRUCTURE FACTORS FOR LARGE cL 



In Section llV Al results from Jimbo and Miwa [8j were used to express the 0(l/cL) correction to the density matrix 
in the form of Toeplitz determinants, which could then be numerically analysed. The method of |8| can also be applied 
to the calculation of the 0(l/cL) correction to the two-point correlation function l|5.1|) . yielding in fact a closed form 
analytic expression. To derive this, we follow the working in and begin by noting that the 0(l/cL) expansion of 
(l2~2|) is 



ip N (xi,x 2 , ■ ■ - ,x N ) 
Now put JV — > JV + 2 (for convenience) and label the particles as 



(l + i- {gHV + » - 1) A - 2 P0 j + . . .) det[e^] 



(6.1) 



< y < xi < . . . < Xj < x < Xj + i < . . . < xn < L. 



(6.2) 



In the definition of gu+2(x,0), x will be the variable as in (|6.2|l and we will take y — > 0. With these labels and 
JV — » JV + 2 the operator in (|6.1|l reads 



-(N + l)^ + ±(-N-l + 2l) i-H-N + l + 2j)l + f (-JV+1+2/) A-2po (± Xl + + y°- 
J i=\ 1=3+1 \i=i " 1 

The determinant in l|6.1[) has the translation invariance property 



(6.3) 



det 



ikjv 



= det 



1 



Jkj(x k -y) 



(6.4) 



since kj = 0. This means that we can write 



dy dx ^-^ dxj 

and using too the fact that we want y — > 0, the operator l|6.3|l reads 



(6.5) 



eta 



\2=1 



' dxi dx 



(6.6) 



In keeping with l|5.1fl . by definition 
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9N+2( 



j(x,0) = lim V* [ dxi . . . dx N \tp N+2 (y,Xi, . . . ,Xj, x,x J+1 , . . . ,x N )\ 2 
where Af 2 is the normalisation, defined by l]2.4[l . and RN,j(y, x) is the region of integration specified by 



(6.7) 



With 



R-N,j{y, x):0<y<xi<...<Xj<x< Xj+i < . . . < a; at < L. 



(6.8) 



(xi,x 2 , ■ • • ,x N+2 ) = (0,xi, . . . ,Xj,x,Xj+i, . . . ,xn) 



(6.9) 



and 



we have 



JV 



A = / , x in — !~ x ^r 

OXi OX 
2 = 1 1 



Jo f) N f) 

2=1 1 i=j+l 



(6.10) 
(6.11) 



lim \i> N+2 (y, x l3 ... , Xj,x, x j+1 ,. . . , x N )\ 2 = | det[e ife ^]| 2 - (| det[e ife ' 2 ']| 2 ) 

y^O c 

+ I 5i (|det[e 1 ^]| 2 )+0^ . 
The normalisation (|2.4() was first expanded in large cL by we display here the first two orders 



(6.12) 



TV 2 = (A^°°)) 2 



l + ___ (A r + i ) + fJ- 

c \cL 



(6.13) 



where (A/^ 00 -*) 2 is the normalisation (|2.4H in the limit cL — > oo. Substituting H6.12|l and l|6.13[) in 1)6. 7fl shows that to 
0(l/cL) 



g N+2 (x, 0) = <^ 2 (x, 0) - ^(JV + l)g^ 2 (x, 0) 



N 

/ dxi...da;ArA(|det[e ife ^ ! ]| 5 

N 

dxi... dx N B,j (| det[e ife ^']| 2 ) 

■ =0 Jr N j{x) 



2po 1_ 

1 1 

c (AT(°°)) 2 ^ 



(6.14) 
(6.15) 

(6.16) 



To proceed further consider | det[e lfej:Ei ]| 2 as a function of xz, (2 = 1, . . . , N). This function vanishes at the three 
points x; = 0, xj, L (j ^ I). 
It follows that 



d 

dxi xi- — | det[ 

R N4 (x) oxi 



; %X(i|2 _ _ 



dxz|det[e ifc ^<]| 2 

R N ,j(x) 



(6.17) 
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while 



d^^-ldet^ 5 ']! 2 = 0. 



R N , 3 ( X ) dx i 



Hence 



g N+2 (x,0) 



ffg&M) ~ ^{N + l)9$ 2 (x,0) - 2 -f (x^ - N^j g { N %(x,0) 



1 1 



c(AA(°°)) 5 



N a r 

V2(j + 1)— / d Xl ...dx N \det[e ik ^}\" 

j=Q ° X J Rn. 3 (x) 



9 ^ 2 (x,0)-^L^ + l)g^ 2 (xAU 



2 1 d 



N 



E- 



(.7 + 1) 



c(^fdx^ o j\(N-j)\J RlfA ^ 



dx\ . . . dxff \ det[e 



where lj6~2*2*|l can be written 



2 1 1 d d ^ - +1 / AT 



dxx... / dxAr|det[e ife ^ ! ]| 2 



2 i i a a 

c (7V(°°)) 2 iV! dx~d£^ 

2i i a ' L 



^fi \[ H f) ^Idet[e ifc ^']| 2 



dxi...da; i v|det[e ifc ^ 1 ]| 2 



c (7V(°°)) 2 AT! 5x i 



we observe that 



1 1 



y dx 1 ...dx A r|det[e ifc ^ 1 ]| 2 = 5 ^ ) 2 (x,0) = p (l-2/(x) 2 ) 



where 



(6.18) 



(6.19) 
(6.20) 
(6.21) 
(6.22) 



(6.23) 
(6.24) 
(6.25) 
(6.26) 



(6.27) 



where we note that (|6.27|l and 
Hence (I6.25|) is equal to 



y(x 



sin 



-iCiV+2,L(x,0) 



(iV + 2)sin(^) po 
also appear in the CUE (Circular Unitary Ensemble), given by Dyson |26j . 



(6.28) 



2 a 



(6.29) 



Regarding H6.26J) . note that 
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_ d_ 



kjxn 1 2 



i / ,l.n detV'< 

x JO J 

N n L 



| [ / ^|det[e ifc ^ ! ]| 2 + ^ / 1) J] / ^ I det[e ife ^] | 

l=l J ° 1=1 J ° = 1,3^1 J ° 



N - x N 

/ dxi TT / dxj|det[e' 
N I dx\ I dx2 ■ 



i=i 



o 



kjXi~\\2 



>0 JO 

Therefore, (|6.26|l is equal to 



dx N \det[e ik ^ l }\ 2 



' dxi [ dxi 



2 1 

c(Af(°°)) 2 {N-l)\dxJ """ X J 



dx N \det[e ik ^ l ]\ 2 



But for the free Fermi system the 3-point correlation function is specified by 



Sa$r+2)( x '°> x d ~ (jV-l)!(jV(°°)) 2 J 



dx 2 ... I dx N \ det[e IKj 
Jo 



ikjxn |2 



1 y(x) y(xi) 
Podetj y( x ) 1 y{x-x x ) 
y(xi) y(x-xi) 1 

Pol 1 - v( x ) 2 - y( x i) 2 - y( x - x if + 2y( x )y( x i)y( x - x i)} 



(6.30) 

(6.31) 

(6.32) 
(6.33) 

(6.34) 

(6.35) 

(6.36) 
(6.37) 



where y(x) is specified by (|6.28JI . Note that this definition is only valid for N > 1, a 3-point correlation function for 
a system of 2 particles does not make physical sense, g^\x, 0, x\) — 0. Now (|6.34|l reduces to 

c^/ dxigff N+2) (x,0,xi). (6.38) 
Finally, adding up all contributions gives the sought closed form expression (here we revert to N + 2 — > N) 



g N (x,0) = g^°\x,0) - ^ 



+ l) g^°\x,0) + - c -^g%°\x,0) + *^J*dx ig <$(x,0,xi) (6-39) 



correct to 0(1/ cL), valid for all N > 2. This has been checked with i|5~9|) and lfS~l3|) . Using (KT27j) and (|6~37l) this 
can be written in the simpler form 



g N (x,0) = ^)(a;,0)+4iV^) { - y(x) [ Po y(x) + y'^)] 



ox 



y(x) / y(x 1 )y[x - x^dx^ 



(6.40) 



which is valid for all N > 3. We have confirmed that this result recovers Ij5.15|l . Now (|6.40(l . and concomitantly its 
structure factor, are readily computed for any N > 3. The structure factor in the limit cL — > oo, is given by 



iV when n = 
S n (N) = \ \ n \/N when < \n\ < N 
1 when |n| > N. 



(6.41) 



17 



In the thermodynamic limit, y(x) becomes 

sinO 7rx) 

y (x) = (6.42) 

p nx 

and now the correlation function follows from 116.40(1 and is (where we use the appropriate scaled variable po / c) 

(sm^ x \ f sm x 

1 j — 4p | _ 3 [ttx cos x + (s — 7r) sinx] 

y{x) J\{x x )y{x-xx)dxi }(^)+o(^) 2 (6.43) 



(9x 



where x = p^irx. This equation was first given by Korepin [27j. Recently in |28j . using the Random Phase Ap- 
proximation (RPA), which they indicated is valid to the po/c correction, the structure factor was calculated in the 
thermodynamic limit from which l|6.43|l is recovered. 

Here we explicitly evaluate the integral in 1)6.43(1 with ((6.4211 

[* sin(x)Si(2x) + cos(z)[Ci(2x) - log(2aQ - 7 ] 

/ y(xi)y{x - xi)dxi — (6.44) 

Jo Ponx 

where 7 is Euler's constant, and Si and Ci are the sine integral and cosine integral functions, respectively. 

In closing this section we note that g^(0,0) in 1(6.40(1 and similarly ((6.431) is zero to leading order, the free Fermi 
result. Deviations from zero do not begin until higher order than 1/cL as can be seen from these equations. 

VII. CONCLUDING REMARKS 

The occupation numbers for the large cL limit given in ((4.21(1 can be continued to the asymptotically large N limit 
giving 

Cn {N)~m+—. (7.1) 

Very strong evidence for the exponent (3 n having the integer value 2 was presented. 

In constructing the ansatz for the occupation numbers, we chose to scale cL by N from the outset. This was done 
because the cL large limit is now compatible with the thermodynamic limit (unlike the small cL limit as discussed 
after 1)8.20(1 ) and further to this we found that the numerical analysis, for any finite N, in constructing the ansatz, 
failed without this scaling. 

The question then is, how does ((7.1(1 compare with the thermodynamic limit. In a very nice work, at a time 
coincident with the seminal work of the Japanese group 0, |29| . the density matrix first given by Lenard in the 
impenetrable limit 0, [31], was extended to the 1 /c correction [30l ]. They used the quantum inverse-scattering 
method in concert with the very important work 29] on the impenetrable Bosc gas in terms of Painleve V theory, 
and obtained 



pM t^f- ( 7 - 2 ) 

X2 *e 



The Fourier transform of 1(7. 2|) gives the momentum distribution, with the Fermi momentum, fcp = Ttpo, 



(7.3) 
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Following the presentation given in [T(j , we can readily see that 1(7.1(1 is one-to-one with 1(7.3(1 . An elementary way 
to see this immediately is to observe that |k| ~ 1/L and thus in the thermodynamic limit |k| ~ 1/iV. Note the integer 
value of 2 for the coefficient of kp / ttc in the exponent 1)7.2(1 and 1)7.3(1 . 

In our work on the impenetrable Bose gas |lOl lllf. we studied the system in all boundary conditions, periodic, 
Dirichlet, Neumann, as well as for the harmonically trapped system. In all cases, we found that in the asymptotically 
large iV limit that ((7.1(1 held in the limit c — > oo. This firmly suggests that the exponent obtained, in this large N 
limit, is universal, being the same for all boundary conditions and (low lying) n modes. Therefore, we anticipate the 
same is true now for our 1(7.1(1 . 

The use of periodic boundary conditions, while facilitating the mathematics, is, nonetheless, a powerful preemptor 
of the analytical properties of the Bose gas, a system of both continuing and immense attraction to the theoretician 
and the experimentalist, in concert. 
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APPENDIX A: BETHE EQUATION SOLUTIONS 
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We have also computed these expansions out to N = 10, due to the inordinate complexities of the numbers, the 
expansions for N > 6 are known only in decimal form. 

The leading term of kj is precisely related to the jth zero of the TVth polynomial, as mentioned in Section [HI Note 
also the universality of the coefficient of the (cL) term, that is —1/24 (as in 1|2.8[) ). 

APPENDIX B: FREDHOLM DETERMINANTS 

Here the Fredholm minor l|4.12|l will be related to a multiple integral, which in turn implies the Toeplitz determinant 
form 1)4.4(1 . Consider the multiple integral 

A N (x,y) = (jf +\J^ \dx 2 ...{j^ +\J^ Jdx N 



N 

x TT 2 sin[(a; — Xj)/L]2 sin[7r(y — Xj)/L] JJ {2 sm[n(xk — Xj)/L]Y 

j=2 2<j<k<N 



(Bl) 



Setting 



g(u) = (l + X X [yl ] )2shx[7T(x-u)/L}2 S m[ ! T(u-y)/L} 



(B2) 



where xf« L = 1 for it e [y, x] and otherwise, allows this to be written 



A N (x,y)= / dx 2 ... dx N Y[g{xi) JJ |e 

, ' jQ 1=2 l<j<k<N 



2iirxk/L 2i-iTXj/L\2 
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By a well known identity (see e.g. Szego [24)) this is equal to a Toeplitz determinant, 



A N {x,y) = (N-l)ldeb 



J j,k=l,...,N-l 

On the other hand the integral in A^(x,y) can be expanded as a power series in A. For this define 



, N {x u x 2l ...,x N )= I] (e^^-e 2 ^^) 

l<j<fc<AT 



and introduce the free Fermi type distribution 
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where Cn,l is such that 



and thus given by 



FF, \ N 

Pn [x,x) = — 



f L f L 
NC N , L = / dxi... dx N \4>(xi,...,x N )Y 
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Now, writing the integrand in l|Bljl in terms of (|B5(I . expanding in A, and making use of the definition 



A™ f x f x 
A N (x, y) = C n ,l dx 2 - ■■ dx n ^ 

n=o n - J v J v 



i+ip ! N F {x,y;x2, ■ ■ ■ ,x n+ i). 

A straightforward calculation (see e.g. Forrester |25|) gives the determinant form 



p N (x,y;x 2 , ■ ■ - ,x n+1 ) = det 



K NL (x,y) i K NL(x,u k )] k=h _ in 



and furthermore shows from (IB8I) that 



C N , L = (N-1)\L N . 
Substituting (|B11(> and (|B10|) in l|B9|) and comparing with (|4.12|) shows 



(N-1)\L N 

With Apf(x,y) in its Toeplitz determinant form i|B4|) . this gives (|4.4|) 



A N (x,0). 
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TABLE III: Values of c ( ( , 0) (iV) and c£ ] (N) for N = 2, 3, 4, 5, 6, 7. Note that this Table extends Table II of 10] 
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TABLE IV: Values of c^ 0) (AT) and c\ L > (AT) for N = 2, 3, 4, 5, 6 
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TABLE V: Values of c { 2 0) (N) and c\ ' (N) for N = 2,3, 4, 5, 6 
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TABLE VI: Values of c^' X) (AT) and c£' 2) (N) for AT = 2, 3, 4, 5, 6 
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TABLE VIII: Values of 4 M) (AT) and (AT) for AT = 2, 3, 4, 5, 6 
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